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ABSTRACT 



We report the discovery of one-hour long tails on the few-minutes long X-ray bursts from the 'clocked burster' GS 1826-24. We 
propose that the tails are due to enduring thermal radiation from the neutron star envelope. The enduring emission can be explained 
by cooling of deeper NS layers which were heated up through inward conduction of heat produced in the thermonuclear shell flash 
responsible for the burst. Similar, though somewhat shorter, tails are seen in bursts from EXO 0748-676 and 4U 1728-34. Only a 
small amount of cooling is detected in all these tails. This is either due to compton up scattering of the tail photons or, more likely, to 
a NS that is already fairly hot due to other, stable, nuclear processes. 
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1. Introduction 

Type-I X-ray bursts, or X-ray bursts in short, result from ther- 
monuclear shell flashes of hydrogen and helium on neutron stars 
(NSs). The fuel is accreted from a Roche-lobe-filling companion 
star. As the accretion of this material progresses, the pressure at 
the bottom of the accreted layer rises to ignition conditions for 
thermonuclear fusion processes like the (hot, /?-decay limited) 
CNO cycle, triple-or process, a-proton capture and the (jfi-decay 
limited) rapid-proton capture (rp) process. In general, the fuel is 
burnt within 1 s in the top few meters of the NS and temperatures 
momentarily reach values as high as a few GK. The layer is cov- 
ered by a non-burning layer, on top of which is the photosphere. 
What one measures is the cooling flux passing through the pho- 
tosphere, with temperatures peaking at about 30 MK. The burst 
duration is primarily determined by the time it takes to cool the 
shell. For H-rich flashes, the duration is further lengthened due 
to prolonged nuclear burning through the rp process. The burst 
duration may range from a few seconds to a few hundred sec- 
onds for the large majority of X-ray bursts. Exceptional dura- 
tions come from very thick helium layers on relatively cold NSs 
in hydrogen-poor ultracompact X-ray binaries (i.e., up to thou- 
sands of seconds; e.g., in 't Zand et al. 2005; Cumming et al. 
2006; in 't Zand et al. 2007) and flashes of very thick carbon 
shells ('superbursts'; e.g., Cornelisse et al. 2000; Cumming & 
Bildsten 2001; Strohmayer & Brown 2002). For reviews on X- 
ray bursts and further references, we refer to Lewin et al. (1993), 
Bildsten (1998) and Strohmayer & Bildsten (2006). 

Sometimes very long tails are seen in X-ray bursts that 
are not related to the aforementioned long bursts. For instance, 
Chenevez et al. (2006) discuss a peculiar burst from GX 3+1. 
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Since the advent of X-ray astronomy, about 100 bursts have been 
detected from this source (e.g., den Hartog et al. 2003) and all are 
shorter than a few tens of seconds except for this peculiar burst. It 
has a prolonged tail that starts at about 25% of the ordinary burst 
peak flux and decays with an e-folding decay time of 1 1 10 s (for 
photons between 3 and 6 keV). Chenevez et al. (2006) hypothe- 
size that the tail is due to rp capture of a rich hydrogen mixture 
that became available after the accretion rate dropped below the 
threshold where hydrogen is burnt in a stable manner (GX 3+1 
was in a 10-yr minimum at about the time of the burst). However, 
the time scales of the slowest /3 decays expected in X-ray bursts 
are at least one order of magnitude shorter than 1110 s (e.g., 
Fisker et al. 2008). A few similar cases (i.e., from systems that 
are clearly not ultracompact X-ray binaries) are described in the 
literature, most notably in Czerny et al. (1987) and Gotthelf & 
Kulkarni (1997). There is a clear duality in the time profile of 
these bursts: they start with an ordinary short-lived burst, fol- 
lowed by a 10 2 ~ 3 s long tail without an unambiguous cooling 
signature. The long tail starts off quite brightly in these bursts, at 
a few tens of percents of the burst peak flux. Since the nature of 
these long tails is not well established, it is worthwhile exploring 
whether there are long tails that start off at a lower fraction of the 
burst peak (i.e., percents instead of tens of percents). 

Burst studies usually concentrate on the brightest parts, 
roughly above 1 % of the peak flux (which is often close to the 
Eddington limit, see Galloway et al. 2008), and for a good rea- 
son. Most bursts come from prolific bursters with mass accretion 
rates above the same 1% level. Since accretion is notoriously 
variable, this makes disentangling burst radiation from accretion 
radiation difficult at fluxes below 1 % of the Eddington limit. 

Exceptions are bursts from NSs that accrete at rates below 
1% of the Eddington limit. This pertains to most persistently ac- 
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creting ultracompact X-ray binaries (UCXBs; a nice example of 
a burst that could be studied with Swift to very deep levels origi- 
nates in A 1246-588; in 't Zand et al. 2008) and bursts from tran- 
sients whose accretion rate has dwindled down to low but non- 
zero values (e.g., in 't Zand et al. 2003). These exceptions show 
tails that are natural extensions of the decays of the bright parts 
of the same X-ray bursts. In other words: there is no prompt/tail 
duality. 

Despite the fact that it is difficult to study bursts at sub- 
1% levels in fast-accreting bursters, the situation is sometimes 
not desperate. This paper presents a study of the unique burster 
GS 1826-24. The accretion is very stable in this source. The vari- 
ability on a time scale of the burst recurrence time (few hours) is 
about 2% rms, so that the recurrence time from burst to burst is 
relatively stable as are the burst peak flux and profile (Ubertini 
et al. 1999; Cocchi et al. 2001; Galloway et al. 2004; Heger 
et al. 2007). This behavior earned it nicknames such as 'the 
clocked burster' (Ubertini et al. 1999) and 'the textbook burster' 
(Bildsten 2000). GS 1826-24 is also notable for a high-energy 
component of the burst emission (i.e., above 30 keV where negi- 
gible amounts of black body emission are expected; in 't Zand 
et al. 1999). Because of its stable bursting behavior and accretion 
rate, GS 1826-24 is excellently suited for low-level burst flux 
studies. The most recent distance determination is 6.07±0.18 kpc 
(for isotropic burst radiation; Heger et al. 2007). 

The existence of long tails in GS 1826-24 was already im- 
plied in in 't Zand et al. (1999) and Thompson et al. (2005). 
The latter work concentrated on the persistent spectrum as mea- 
sured with Chandra and RXTE. This emission could be success- 
fully modeled by a combination of two comptonization compo- 
nents (see also Thompson et al. 2008) due to the presence of hot 
plasmas in the immediate neighborhood of the NS. Bursts were 
studied as well, and their spectra could be modeled by a black 
body component plus the variation of one of the persistent comp- 
tonization components. This is in line with the high-energy burst 
emission seen by in 't Zand et al. (1999) and in contrast with the 
so-called 'standard' modeling of burst spectra where the burst 
emission is modeled solely by one black body component and 
the persistent spectrum is unaffected by the burst emission. A 
simple physical interpretation of the change of one comptoniza- 
tion component would be that a hot plasma up-scatters some of 
the burst thermal photons to higher energies and itself is cooled 
down by the soft photons. This model was applied to 1000 s of 
burst data in Thompson et al. (2005). The final 850 s of the data 
were modeled by comptonization only. 

In this paper we make a study of the long tail in GS 1826- 
24 employing all RXTE data available and cross checking with 
XMM-Newton data. Thus, we obtain a superior statistical qual- 
ity and are able to probe the tail longer than Thompson et al. 
(2005). We supplement this analysis with briefer investigations 
of two other prolific bursters, EXO 0748-676 and 4U 1728-34. 
We propose an explanation for the tail. 

2. Observations 

We use data from the Proportional Counter Array (PCA) on the 
Rossi X-ray Timing Explorer (RXTE) and the European Photon 
Imaging p-n junction camera (EPIC-pn) on the XMM-Newton 
observatory. These data sets were chosen because they are com- 
plementary in photon energy bandpass (2-60 and 0.1-12 keV, 
respectively) with similar sensitivities and because they detected 
more than 10 X-ray bursts each. 

The PCA comprises 5 identical proportional counter units 
(PCUs) with a total net collecting detector area of about 



Table 1. Selection of X-ray bursts from GS 1826-24 studied here. For 
more details, see Galloway et al. (2008) 



MJD 


ObsID 


Times 


Active 


A" 






covered wrt 


PCUs 


(') 






burst start 










time (sec)* 






50971.23019 


30054-04-02-01 


-2600/+1200 


0,1,2,3,4 


0.0 


50971.70032 


30054-04-02-000 


-2700/+900 


0,1,2,3,4 


0.0 


50988.82503 


30060-03-01-01 


-1700/+1800 


0,1,2,3 


0.0 


51725.87648 


50035-01-02-00 


-3000/+200 


0,2,3,4 


0.5 


51 726.7 1988 c -"-' 


50035-01-02-02 


-1200/+2100 


0,1,2,3 


0.5 


51728.77338'-"-' 


50035-01-02-04 


-4000/+2600 


0,1,2,3 


0.5 


51814.34546 c - w 


50035-01-03-10 


-2200/+1500 


0,2 


0.5 


52484.41775 c - w 


70044-01-01-000 


-4000/+2200 


0,2,3 


0.0 


52484.56684 c 


70044-01-01-00 


-2400/+700 


0,2,3 


0.0 


52485.00672^-' 


70044-01-01-02 


-4000/+2600 


0,2,3 


0.0 


52537.22347 


70025-01-01-00 


-4000/+2400 


0,1,2,3 


0.0 


52738.47686 c - w 


80048-01-01-04 


-4000/+3200 


0,2,3 


0.0 


52738.74636 c -"-' 


80048-01-01-07 


-4000/+2700 


0,1,2,3 


0.0 


52820.56256 c - p -' 


80049-01-01-00 


-1100/+2000 


0,2,3,4 


0.0 


52893.05439 


80105-11-01-00 


-4000/+3000 


0,2,3 


0.0 


53205. 19659 c - w 


70044-01-02-00 


-1900/+1500 


0,2,3 


0.0 


53207.07579 


90043-01-01-01 


-4000/+3600 


0,2 


2.5 


53207.2191 l c 


90043-01-01-01 


-4000/+1400 


0,2,3 


2.5 


53206.49150 c - w 


90043-01-01-020 


-4000/+2800 


0,2,3 


2.5 


53206.63615 CJV 


90043-01-01-02 


-4000/+1900 


0,2,3 


2.5 


53956.08413"-' 


92031-01-01-000 


-4000/+2200 


0,2,4 


13.8 


53956.22695"-' 


92031-01-01-000 


-2100/+1100 


0,1,2 


13.8 


53957.06785"-' 


92031-01-01-010 


-4000/+2000 


0,2,3 


13.8 


53957.20458"-' 


92031-01-01-010 


-1800/+1600 


0,2,3 


13.8 


53959.15478"-' 


92031-01-02-000 


-4000/+2900 


1,2 


13.8 


54167.72354 c -"-' 


91017-01-01-000 


-1900/+1100 


1,2 


0.0 


54168.31443 c -"-' 


91017-01-02-03 


-2000/+1600 


2 


0.0 


54168.75518 CJV 


91017-01-02-010 


-4000/+2500 


0,2,3 


0.0 


54168.90229 c -"-' 


91017-01-02-010 


-1900/+1600 


0,2 


0.0 



" Off-axis angle; 6 These values are rounded off to an integer times 100 
s and may include gaps. c The spectra for these ObsIDs were used to 
generate average time-resolved spectra (cf., Fig. 4); "Data were used 
to calculate average pre-burst power density spectrum (Fig. 6); 'Data 
were used to calculate average tail power density spectrum. 



8000 cm 2 in the 2-60 keV band (3-20 keV is well calibrated; 
Jahoda et al. 2006). During any observation any number of PCUs 
between 1 and 5 are on. The PCA is a non-imaging device 
with a spectral resolution of 20% (full width at half maximum; 
FWHM) and a 2 degree wide field of view (full width to zero 
response). There are no other bright X-ray sources in the field 
around GS 1826-24. 

RXTE is particularly suited for the study of low-mass X-ray 
binaries, and since its 1995 launch has accumulated unprece- 
dented exposure times on many X-ray bursters, one of them 
being GS 1826-24. Galloway et al. (2008) compiled a cata- 
log of all bursts detected between Feb 8, 1996, and June 3, 
2007. The total PCA exposure time on GS 1826-24 over 127 
observations in this time frame is 929 ks. 65 bursts were de- 
tected between Nov 5, 1997, and March 10, 2007, and the mean 
burst rate was 0.25 hr _1 . None of the bursts exhibit evidence 
for photospheric radius expansion, so that the peak luminos- 
ity must have been below the Eddington limit. The bolomet- 
ric absorption-corrected peak flux ranged between 23.5 + 0.8 
and (29.6 ± 0.8) x 10~ 9 erg cm _2 s _1 and the fluence between 
0.923 + 0.016 and (1.059 + 0.004) x 10" 6 erg crrr 2 (Galloway 
et al. 2008). The decay can be modeled by two subsequent expo- 
nential decays. The e-folding decay times ranged between 12.9 
and 23.3 s for the first decay and between 40.5 and 57.2 s for 
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Fig. 1. Average time profile of 29 X-ray bursts from GS 1 826-24 as 
measured with PCU2 on RXTE/PCA, in two photon-energy bandpasses 
and at logarithmically scaled resolution. 



the second decay. The average time scale for the decay (defined 
as the fluence divided by the peak flux) is between 30 and 45 s, 
which implies that the 0.1% level is reached within 400 s. 

XMM-Newton observed GS 1826-24 on two occasions: 
starting on April 6, 2003, for 108 ks and April 8, 2003, for 92 ks. 
A complete account of these observations is provided in Kong 
et al. (2007). Nine X-ray bursts were detected in the first and 
seven in the second observation. The final bursts in each obser- 
vation suffered from high background levels and were excluded 
from the analysis. All X-ray detectors were on, but we concen- 
trate on the 0.1-12 keV EPIC-pn measurements because that in- 
strument (Striider et al. 2001) is by far the most sensitive for 
our analysis. It has an effective area that ranges between 1000- 
3000 cm 2 for 0.5-2 keV and 900 cm 2 at 2-6 keV. The spectral res- 
olution is 2.5% (FWHM) at 6 keV. The instrument was switched 
to Fast Timing mode, implying that the central CCD (of the 12 
available), encompassing 13 '6x4 .'4 of the field of view, is read 
out 1-dimensionally (along the 4.'4 side) at a 1.5 ms resolution. 
Source photons were extracted between RAWX values of 30 and 
45, with pixel patterns below 5 and grade 0; background pho- 
tons were extracted between RAWX=10 and 25. We refer to Kong 
et al. (2007) for further details, noting that SAS version 7.1.0 
was employed for our data analysis. 

3. Data analysis 

We selected 29 out of the 65 PCA-detected bursts that have data 
available between 4000 and 1000 s before the burst start and 
between 400 and 3000 s afterwards, see Table 1 . The pre-burst 
interval was based on the general trend that the flux was lowest 
there and should provide the best estimate of the accretion flux. 
We processed the data of the 29 bursts to an average light curve 
as follows. Taking the 'standard-products' light curves as start- 
ing point (these are data from PCU2 in 5 energy bands, corrected 
for particle-induced background and with a resolution of 16 s) 
we determined the pre-burst flux level from data between 4000 
and 1000 s prior to the burst (note that data gaps are common in 
this 3000 s time frame), subtracted that from all flux measure- 
ments, co-aligned the data at the start time of the burst (defined 
as the time when the photon flux exceeds 25% of the peak flux; 
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GS 1826-24 / RXTE-PCA / 2-9 KEV / 29 BURSTS 
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Fig. 2. Average time profile of 29 X-ray bursts from GS 1826-24 as 
measured with PCU2 on RXTE/PCA in 2-9 keV and at logarithmically 
scaled resolution, together with an exponential decay function (smooth 
curve) as fitted between 600 and 5500 s (r = 1252 s; see Table 2) 



GS 1826-24 / XMM-EPIC/PN / 2-9 KEV / 14 BURSTS 
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Fig. 3. Average time profile of 14 X-ray bursts from GS 1826-24 as 
measured with XMM/EPIC, in one bandpass and at a logarithmically 
scaled resolution. The count rate scaling is different from Fig. 1 due to 
the different time resolution at the peak and the occurrence of saturation 
effects in the XMM-Newton data. The smooth curve shows an exponen- 
tial function, fitted between 600 and 3000 s. The decay time was found 
to be t = 1006 ± 26 s. The goodness-of-fit is formally unacceptable 
(X 2 /v = 179/22), so this result is only indicative. 



Galloway et al. 2008) and averaged all bursts. Figure 1 shows 
the resulting profiles in 2 to 4 and 4 to 9 keV, normalized to the 
peak flux. Apart from the well-known burst profile with a dura- 
tion of about 400 s, it shows the striking appearance of an addi- 
tional burst component that lasts approximately ten times longer 
at flux levels between 10~ 3 and 10~ 2 times the burst peak value. 
A comparison with higher time-resolution data shows that the 
averaging of the peak in 16 s time bins lowers the peak flux by 
approximately a factor of 2. While the initial decay of the profile 
shows the classical cooling of the bursts (the high-energy flux 
decaying faster than the low-energy flux), the slow decay does 
not show obvious cooling. 
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Fig. 4. Modeling results for time-resolved PCA burst spectra of aver- 
age of 17 bursts. The model consists of a black body and power-law 
component absorbed by a fixed column of Nu = 3.1 x 10 2 ' cirr 2 . The 
panels show from top to bottom: the bolometric flux for the black body 
component (solid curve) and the unabsorbed 3-20 keV flux found for 
the power-law component (dashed curve; gaps indicate times when the 
power law is not detected), the black body temperature, the emission 
area in terms of the radius of a sphere at a canonical distance of 10 kpc, 
and the goodness-of-fit in terms of reduced x 2 (the number of degrees 
of freedom is 17 or 18). The black body temperature is, beyond 150 s, 
fixed to the values found in the XMM-Newton data. 



The data can be described satisfactorily by an exponential 
decay function between 600 and 5500 s. For the 2-9 keV time 
profile the fit is shown in Fig. 2. the decay time is 1252 + 25 s 
(X 2 /v = 16.9/13). Resolved in the two bands the e-folding decay 
times are 1261+29 s in 2-4 keV (x z /v = 16.3/10) and 1381+29 
in 4-9 keV (as measured between 600 and 3500 s after burst 
onset for lack of statistics beyond 3500 s; x 2 /v = 53.4/10). Since 
it is expected that cooling by thermal/photon diffusion follows 
a power law (Eichler & Cheng 1989a) we fitted such a law and 
find, between 300 and 2300 s, a power law index of -0.93 +0.02. 

We sought verification of the long tail in data from the 
XMM-Newton observations. Figure 3 shows the 2-9 keV light 
curve, averaged over 14 out of the 16 X-ray bursts (leaving out 
the last burst of each observation, being compromised by in- 
creased background levels). The same long tail is seen as with 



Fig. 5. Modeling results for time-resolved 0.7-10 keV burst spectra of 
14 bursts detected XMM-Newton observation. Spectral channels were 
binned so that each bin contains at least 15 photons. The number of 
degrees of freedom ranges between roughly 500 and 2000 (for 2 spectra, 
one for each observation) which explains the small scatter in xl- The 
anomalous values at 20 s are due to data drop outs due to full data 
buffers from high photon rates (see Kong et al. 2007). 



RXTE. The e-folding decay time is t = 1006 + 26 s between 
600 and 2400 s which is similar as for the RXTE curve. A tail is 
also seen in data below 2 keV, but the statistics are not so good 
to reveal it beyond 10 3 s. The e-folding decay time in this low- 
energy bandpass is 524 + 172 s between 400 and 1000 s, which 
at least shows that it is not longer than for the 2-9 keV band. 

We performed time-resolved spectroscopy. Spectra were ac- 
cumulated for 17 bursts that have small off-axis angles in the 
PCA and have identical detector-voltage settings (i.e., they are 
said to be in the same 'gain epoch'; see specification in Table 1). 
Only PCU2 data were employed to obtain a homogeneous data 
set. The time resolution of the spectral extraction was chosen 
to vary between 1 s early on in the burst and 200 s at the end. 
Each burst was divided in 72 time bins that are identical with 
respect to the burst start times. Spectra were extracted from 
event mode data with s/w tool seextrct for these time intervals 
and corrected for the particle-induced background as determined 
with pcabackest in the same time interval. For each burst, a 
pre-burst spectrum was generated from data available between 
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Fig. 6. Geometrically binned Fourier power density difference spectrum 
between pre-burst and tail data (in formula: tail pds minus pre-burst pds 
divided by pre-burst pds). 

4000 and 1000 s prior to the burst, which also was corrected for 
particle-induced background. Subsequently the spectra for all 17 
bursts were averaged in their respective time frames, and the av- 
erage pre-burst spectrum was subtracted from the 72 burst spec- 
tra. These spectra were, between 4 and 20 keV, modeled with 
a simple absorbed black body function, employing a fixed ab- 
sorption column of 7V H = 3.1 x 10 21 ctrT 2 (see XMM-Newton 
analysis below). The photo-absorption cross sections were taken 
from Balucinska-Church & McCammon (1992) and the compo- 
sition of the absorbing material from Wilms et al. (2000). About 
half of the spectra (in the brightest phase) turn out not to be con- 
sistent with this model and an additional power-law component 
with a fixed index (equal to that of the pre-burst data) was in- 
cluded in the model. The fitted values for the various black body 
parameters are shown in Fig. 4. All time intervals are well fit- 
ted. During the first 100 s the black body parameters are as ex- 
pected for an X-ray burst: a temperature peaking at an equivalent 
of 2.3 keV and gradually decreasing after that, and an emission 
area that remains approximately constant after the rise phase. 
However, after 100 s the picture changes: the inferred emission 
area decreases sharply by at least an order of magnitude. The 
bolometric unabsorbed fluence in the 300-1500 s time frame is 
3.4xl0~ 8 erg ctrT 2 . This is about 3% of the fluence in the prompt 
burst (Galloway et al. 2008). 

A similar time-resolved spectral analysis was performed 
XMM-Newton data from the 14 low-background bursts, except 
that this involves an analysis of 0.7-10 keV photons and chan- 
nels were binned to make sure that the number of photons per 
bin was in excess of 15 to ensure applicability of the x 1 statis- 
tic (such a procedure was not necessary for the RXTE data). 
Prior to the burst spectral modeling we modeled the pre-burst 
data to find Nh and the photon index to apply to the burst data. 
These are N H = 3.1 x 10 21 ctrT 2 and F = 1.47 (jc 2 v = 1.001 for 
v = 2923 over two spectra, for a systematic uncertainty per bin 
of 2%). The results of the burst spectra modeling (Fig. 5) are 
generally consistent with the RXTE results, except at the start 
and at the end of tail (beyond about 1000 s after burst onset) 
where a comparison becomes difficult because of statistical is- 
sues. This shows that the drop in radius at 100 s as seen with 
RXTE is not related to the lack of low-energy coverage of that 
instrument since it is also seen with XMM-Newton for which the 
bandpass is extended with the 0.7-4.0 keV photon energy range. 



10.000 r 




Time [sec] 



Fig. 7. Modeling results for time-resolved PCA burst spectra of average 
of 36 bursts from 4U 1728-34. The model consists of a black body 
and power-law component absorbed by a fixed column of Af H = 2.6 x 
10 22 cirT 2 . The panels show from top to bottom: the bolometric flux for 
the black body component, the black body temperature, the emission 
area in terms of the radius of a sphere at a canonical distance of 10 
kpc, and the goodness-of-fit. The data points beyond 250 s have large 
X 2 values (outside range of bottom plot) and do not fit the model. This 
is most likely due to variability in accretion flux. 

We studied the timing properties of the tail in comparison to 
those of the pre-burst data. Fourier power density spectra were 
generated of data in the pre-burst -4000/- 1000 s time frame and 
in the tail +1000/+3000 s time frame. Event mode data were em- 
ployed at a time resolution of 2~ n s (roughly 0.5 ms), as long as 
they were available (see Table 1). Power spectra were made in 
1968 8-s data stretches for the pre-burst data and averaged, and 
for 2574 8-s data stretches in the tail data and averaged. Out of 
the FTOOL package vs 6.5, powspec was employed for this pur- 
pose. The relative difference between these two average power 
spectra is presented in Fig. 6. There is no difference between 
both spectra, indicating that nothing significantly changed in the 
accretion stream between these periods. 

4. Other bursters 

We checked the average burst profiles of other known prolific 
bursters in the PCA data archive, employing the RXTE bursts 
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catalog (Galloway et al. 2008). Most of these have variable ac- 
cretion fluxes and are, therefore, difficult to analyze for the rea- 
sons discussed in the introduction. Nevertheless, they sometimes 
indicate long tails, although never as long as in GS 1826-24. In 
this section we report briefly two cases. 

4.1. 4U 1 728-34 = GX 354-0 

4U 1728-34 was observed 346 times for a total PCA exposure 
time of 1.94 Ms and a total of 106 bursts were detected; a large 
portion of these, 69, show photospheric radius expansion. All 
bursts are short and have time scales between 4.4 and 8.7 s 
(Galloway et al. 2008). 40 bursts have good coverage and are 
not within too wild accretion flux variations to search for a tail. It 
turns out that the average profile of these bursts extends to about 
800 s which is 10 2 times longer than the initial burst phase. The 
exponential decay times are 477+80 s for 2-4 keV, 314+26 s for 
4-9 keV (between 300 and 800 s) and 306+12 s for 2-9 keV. A 
power-law fit to the tail in the latter bandpass between 100 and 
800 s yields a decay index of 1.17 ± 0.02. 

We performed time-resolved spectroscopy in a similar man- 
ner as for GS 1826-24 (i.e., the persistent flux as determined 
from pre-burst data was subtracted prior to the spectral mod- 
elling). There was no need to include a power law. 36 bursts 
were selected for this procedure and only PCU2 data were em- 
ployed. Since there is no majority within one RXTE PCA gain 
epoch, we employed bursts from epochs 3a, 3b, 4 and 5, aver- 
aged them per gain epoch, modeled per gain epoch and averaged 
the fitted parameters over the epochs per time bin. A warning 
is appropriate here: the bursts in 4U 1728-34 vary more than in 
GS 1826-24; averaging them will smooth out short features. The 
results yield the profiles in Fig. 7 (xl was averaged as well; al- 
though that number does not adhere to x 2 statistical properties, 
it does give a sense of the overall quality of the fit per bin). The 
bolometric flux can be followed downward over 3 decades. It 
has a smoother evolution than for GS 1826-24, possibly because 
the Eddington limit is reached for most bursts so that the flux 
flattens against a ceiling representing that limit. Furthermore no- 
ticeable is the drop in radius beyond 10 s, a similar effect as seen 
in GS 1826-24 beyond 100 s. The fluence of the tail is 9% of that 
of the prompt emission, taking 100 s as the boundary between 
prompt burst and burst tail. We scale the boundary with the de- 
cay time T2 of the latter part of the prompt burst (see Table 2, as 
read from Galloway et al. 2008). 

4.2. EXO 0748-676 

EXO 0748-676 is well known for the accretion-disk edge-on line 
of sight, causing eclipses and dips every orbital period of 3.8 
hr (Parmar et al. 1986), and also for exhibiting very short burst 
recurrence times in certain accretion rate regimes. Boirin et al. 
(2007) discovered with XMM-Newton that there are times when 
this source exhibits 3 bursts in a row within only half an hour. 
The first burst in such a 'triplet' always shows a longer tail than 
the subsequent bursts. In the average burst profiles the e-folding 
decay time is about 2.5 times longer for the first bursts than for 
the subsequent bursts (50-55 s versus 14-19 s). 

RXTE observed EXO 0748-676 ninety-four times for a 
total exposure of 1.39 Ms. Eighty-four bursts were detected. 
However, many of these observations were concentrated on 
catching eclipses and, thus, involve only small time stretches 
rendering burst tail studies impossible. Also, the bursts are rather 
weak so that multiple PCUs are needed to perform meaningful 



EXO 0748-676 / PCUs 0+2+3 / 1 6 s bins / Channels 1 -1 3 




-1000 -500 500 1000 1500 



Time [s] 

Fig. 8. Light curve of triple burst from EXO 0748-676 seen on 
September 12, 2006, at 16:30:36 UT (ObsID 92013-01-03-000). The 
peaks of the bursts are truncated to zoom in on the long tail of the first 
burst. The peak fluxes of the 3 bursts are 785.5, 462.1 and 86.8 c/s. The 
horizontal line is plotted as a reference for the out-of-burst flux. 



analyses. Rather, we employ XMM-Newton data to determine 
average burst profiles. EXO 0748-676 is the burster which, with 
XMM-Newton, was most intensely observed. In 158 hr of ex- 
posure investigated by Boirin et al. (2007) there are 33 singlet 
bursts detected, next to 14 doublets and 5 triplets. We deter- 
mined average time profiles of 3 1 singlets that are not affected 
by eclipses or dips within 2000 s, in the same bandpasses as 
for the RXTE data for GS 1826-24 and 4U 1728-34. These also 
show a long tail. The exponential decay times are 266 + 122 s 
and 260 ± 120 s for 2-4 and 4-9 keV respectively. Again no clear 
cooling is observed in this tail. The equivalent power-law decay 
index is 2. 1 ± 0.8 for 2-4 keV and 1 .6 + 0.8 for 4-9 keV. 

RXTE data do show at least one interesting burst. It is a burst 
triplet that occurred on September 12, 2006 (onset of first burst at 
16:30:36 UT). The light curve, as extracted from standard 2 data 
in ObsID 92013-01-03-000, is presented in Fig. 8. Standard-2 
data from PCUs 0, 2 and 3 were added within 2.0-7.3 keV and 
the time resolution of these data is 16 s. The results of time- 
resolved spectroscopy are shown in Fig. 9. The unabsorbed bolo- 
metric fluence ratio over the 3 bursts is 5 : 1 : 1/10 (for the first 
burst accumulating fluence over only the first 300 s, so excluding 
the long tail). The light curve clearly shows a long tail to the first 
burst. The e-folding decay time is 300 ± 42 s (measured between 
400 and 1200 s after burst onset and excluding the prompt flux 
of the 2nd burst). This is similar as found in the average XMM- 
Newton burst profile. The fluence in the tail, excluding the three 
bursts, is 4.5xl0~ 9 erg crrT 2 which is 1.2 + 0.1% of the prompt 
emission from the first burst (taking 300 s as the boundary be- 
tween both). It is interesting to note that, in the bolometric flux, 
there is no clear boundary between the long tail in this triple 
burst and the prompt burst, in the sense that there is no excess of 
the early emission above the backward extrapolation of the tail. 
This is contrast to the situation for GS 1826-24 (cf., Fig. 4 top 
panel). The same applies to 4U 1728-34 (Fig. 7). 

What is most remarkable about this long tail is that it is seem- 
ingly unaffected by the occurrence of the second burst. The tail 
progresses undisturbed along the same decay curve. This sug- 
gests that the first burst and its tail emission have a different ori- 
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Fig. 9. Time-resolved spectroscopy of the triple burst from EXO 0748- 
676. A pre-burst spectrum was subtracted before fitting the spectra with 
an absorbed black body model with N H = 8 x 10 21 cirr 2 . For data 
between 60 to 1200 s, we fixed kr for the tail emission to 1.4 keV, equal 
to the average value if kr is left free for the relevant time intervals. 



gin in the NS envelope than the second burst, either a different 
layer or a different locality on the surface. However, this is not 
testable: the fiuence of the second burst is five times smaller than 
that of the first. If the tail would be proportionally smaller, the 
statistical significance of it will have dropped below the detec- 
tion threshold. 



5. Discussion 

Our measurements are summarized in Table 2. We find that X- 
ray bursts from GS 1826-24 show a long tail. In other words, 
they exhibit a dual time profile with a prompt burst phase lasting 
a few hundred seconds and a tail phase with an e-folding decay 
time of 10 3 s, a flux level of less than «1% of the peak flux and 
a fluence that is 3% of that of prompt burst. Beyond the first 
200 s the spectrum can be modeled by a black body with a very 
slowly decreasing temperature of «0.9 keV and a strongly de- 
creasing emission area. We find similar tails in RXTE data of 4U 
1728-34 and RXTE and XMM-Newton data of EXO 0748-676, 
although the contrast between prompt and tail phase is less pro- 
nounced in those cases. Generally the detection of tails is diffi- 
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Fig. 10. Average observed (red histogram) and predicted (solid) time 
profiles of bolometric flux, after normalization at the respective peak 
values and alignment at the peak. The vertical dashed lines indicate 
the times for the depth profiles plotted in Fig. 11. The inset shows the 
model on a broader and linear time scale, exhibiting the slow luminosity 
increase between bursts from burning of increasing amounts of accreted 
hydrogen in the hot CNO cycle. 



cult because of confusion with variable accretion radiation. Still, 
bright pronounced tails have been reported in the literature for a 
few individual bursts. The two most obvious questions about the 
long tails are: what is the physical cause and why is hardly any 
cooling detected? We investigate the first question in Sect. 5.1, 
keeping in mind only GS 1826-24. In Sect. 5.2 we touch on the 
second question, considering 4U 1728-34 and EXO 0748-676 as 
well. In Sect. 5.3 we compare GS 1826-24 with the other cases 
listed in Table 2. 

5. 1 . Origin of the long tail 

What is the origin of the 10 3 s time scale of the long tail? One 
idea spawns from the detection of a hot plasma surrounding 
the NS (Thompson et al. 2005), namely that prompt burst pho- 
tons are trapped in the hot plasma through scattering and that 
the 10 3 s time scale is the time it takes to drain the plasma of 
those photons. The maximum size of the plasma in a 2. 1 hr orbit 
(Homer et al. 1998) is of order 10 11 cm. For an optical depth of 6 
(Thompson et al. 2005), this implies a drainage time of at most 
10-10 2 s. This is one to two orders of magnitude shorter than 
observed. Therefore, if 2.1 hr is indeed the orbital period (this 
needs to be verified; e.g., Mescheryakov et al. 2004), the 10 3 s 
time scale cannot be explained by this idea. 

Another idea is that the tail is due to cooling of layers that 
are deeper than the flash layer. In a thermonuclear shell flash, 
heat is transported upward (to be radiated by the photosphere) as 
well as downward (heating up deeper layers through conduction; 
Eichler & Cheng 1989b; Strohmayer & Brown 2002). Deeper 
layers, extending down to the crust, consist of the ashes of the 
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Table 2. Summary of measurements 



Burst source 


Mean burst duration (s) fl 


Tail/prompt ratio 


Ref." 




time scale c 


d 
'?. 


tail* 


in time^ 


in fiuence 8 






Source^ 


discussed here 








GS 1826-24 


39 .4 ± 0.2 


44.6 


1252 ± 25 


32/28 


0.03 (300 s) 










(600-5500 s) 














^ 2 /v=16.9/15 








4U 1728-34 


5.85 ± 0.02 


12.6 


313 ± 17 


53/25 


0.09 (100 s) 










(200-800 s) 














X 2 /v= 13.4/11 








EXO 0748-676 


13.7 ±0.1 


36.3 


300 ± 42* 


22/8 


0.01 (300 s) 










(400-1200 s) 














^ 2 /v=141/44 








Literature cases 


GX3+1 


3 + 1 


6 


1110±170 


370/185 


40 (10 s) 


1 


M28 




7.5 


800-3250 


~ 100 


> 1.4' (90 s) 


2 


AqlX-1 


18.4 




627 + 100 


34 




3 


OSO-8 long burst 


44.8 




215.7 


5 


0.01 (150 s) 


4 



"Averaged over the bursts cataloged by Galloway et al. (2008); ^References: 1 - Chenevez et al. (2006), 2 - Gotthelf & Kulkarni (1997), 3 - 
Czerny et al. (1987), 4 - Swank et al. (1977); The time scale is defined as the burst fiuence divided by the peak flux; d T 2 is the e-folding decay 
time of the final part of the prompt burst; The time range in parentheses refers to the interval in which the fit was carried out; The two numbers 
refer to division by the mean time scale and t 2 , respectively; s Fluences for GS 1826-24 and 4U 1728-34 were determined from data averaged 
over multiple bursts; that for EXO 0748-676 from the RXTE triple burst. The time between parentheses refers to the interval for the prompt 
emission; ''this value applies to the triple burst only and excluding times for the second burst ; 'based on counts of 0.8-12 keV photons in 

ASCA-GIS (table 1 in Gotthelf & Kulkarni 1997) 



nuclear H and He burning and are rich in elements up to a mass 
number between 60 and 100 (Schatz et al. 2001; Fisker et al. 
2008). Radiative cooling of these deeper layers can explain the 
long tail. The time scale of the long tail points to a column depth 
that is 10 to 30 times larger than that of the burning layer. The 
amount of fiuence in the tail is one to two orders of magnitude 
smaller than in the prompt burst, suggesting less heating down- 
ward than upward. 

The idea that the long tail results from cooling of deeper lay- 
ers is corroborated by model calculations. Heger et al. (2007) 
calculated various sequences of flashes specifically for GS 1826- 
24 with different mass accretion rates and metallicities. Their 
model 'A3' is the one whose recurrence time of 3.85 hr matches 
best the majority of our bursts. This model assumes a mass ac- 
cretion rate of 1 .58x 10~ 9 M G yr~' (or 1 x 10 17 g s _1 ) and a metal- 
licity of Z = 0.02. The long-term time profile of the radiated 
luminosity is very similar to Fig. 27 of Woosley et al. (2004). 
Ignoring the period before the first burst and after the accretion 
turn-on, the interburst time profile is characterized by a grad- 
ual decline till about 3000 s after the burst onset followed by a 
gradual increase for about 10 4 s untill the next burst. The decline 
is due to cooling of the deeper layers; the increase is due to hot 
CNO burning of newly accreted hydrogen in the burning zone. In 
Fig. 10 we plot the observed bolometric flux (see also Fig. 4) and 
the average luminosity profile of 29 bursts from the 'A3' model 
by Heger et al. (2007). The fluxes and luminosities were normal- 
ized to the peak value, and the light curves were aligned at the 
peak. Times in the model were corrected for general relativistic 
effects (through a multiplication with 1 +z = 1.26; Woosley et al. 
2004). The time profiles are an excellent match, all the way to 
the long tail. Figure 1 1 shows for this model the evolution of the 
depth profiles of temperature, net outward luminosity and nu- 
clear energy generation rate per unit mass. Figure 12 presents 
a more detailed view of the latter panel, showing the dynamic 
depth profile of the specific nuclear energy generation rate, an- 
notated with the various nuclear processes playing a dominant 



role at the various locations. Most of the nuclear burning occurs 
for column depths y < 2 x 10 8 g cm 2 , but layers are heated that 
are 10 times deeper (see first panel of Fig. 1 1). The energetically 
less important nuclear burning below 2 x 10 8 g cirT 2 ('heated 
(a, yY or heat-induced a capture by heavier isotopes) is actually 
the result of conductive heating from the shallower layers. 

The inward heating is due to conduction. To first order, 
the heat transport scales with A/Z 2 (Yakovlev & Urpin 1980; 
Bildsten & Cutler 1995; Cumming & Bildsten 2001) where Z 
is nuclear charge and A nuclear mass number. Therefore, one 
may expect less inward heating for a layer with heavier iso- 
topes. The composition depends, on its turn, on the composition 
of the donor atmosphere and the accretion rate. Thus, one may 
expect more inward heating in for instance ultra-compact X-ray 
binaries or for accretion rates and burst regimes with high a in 
which most hydrogen is burnt through the hot CNO cycle instead 
of the rp process, because the rp process produces the heavi- 
est elements. Unfortunately, the A/Z 2 proportionality of the heat 
transport is only a crude approximation so that inferences are not 
straightforward, for instance by comparing different bursters or 
bursts from the same source at different accretion rates. 

There are other dependencies of the conductivity as well, 
such as on ignition depth which may vary from source to source 
and burst to burst. This will not only have an effect on the dura- 
tion of the tail, but also on the fiuence ratio between the tail and 
the prompt burst. For a relatively shallow ignition, the inward 
heating will not go as deep and the tail will be short and less 
fluent. This may explain qualitatively why only the first burst in 
the triplet from EXO 0748-676 has a long tail that is unaffected 
by subsequent bursts. The subsequent bursts ignite at shallower 
depths. 

As mentioned above, the 'A3' model predicts a trend in the 
NS luminosity between bursts consisting of a gradual decline for 
3000 s after a burst followed by a gradual increase for 10 4 s up 
to the next burst (see Fig. 10). In practice it is difficult to dis- 
entangle the NS luminosity from the flux measurements, but it 
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Fig. 11. Profiles of temperature (top panel), net outward luminosity (middle; dashed parts indicate inward luminosity) and specific nuclear energy 
generation rate (bottom) calculated as a function of depth for a flash model specifically for GS 1826-24 (model 'A3' of Heger et al. 2007). Seven 
profiles are shown. The black curve is for t = (burst onset), the other curves are for the times indicated with vertical lines in Fig. 10 with the same 
color. The top panel shows the inward heating (compare with the location of the heat source in the bottom panel, down to only 9 x 10 8 g cm -2 ) 
that is thought to be responsible for the long tails on the X-ray bursts of GS 1826-24. See also Fig. 12. 



interesting to show the straightforward observed photon rate, see 
Fig. 13. This is the same kind of measurement as shown in Fig. 2, 
except that it is shown linearly and between 12,000 s before and 
after the burst time. To obtain as much data as possible, particu- 
larly far away from the bursts when the data coverage is less than 
close to the bursts, we included 27 additional bursts for a total 
of 56. This plot does not show a clear increasing trend except 
for the final 2xl0 3 s before the next burst. A similar behavior, 
but of worse statistical quality, is apparent from XMM-Newton 
data. One would expect a linearly increasing trend for the hot 
CNO cycle according to L = 4 x 10 34 (Zc N o/0.02)Mi 7 f hr erg s _1 , 
with Zcno the CNO abundance, M\i the mass accretion rate in 
10 17 g s and V time in hours, because the amount of accu- 
mulated fuel grows linearly with time and the hot CNO burning 
rate is a constant that depends only on the CNO mass fraction 
(Hoyle & Fowler 1965). The maximum slope consistent with 
the data shown in Fig. 13 is, for M\i = 1 (Heger et al. 2007) and 
assuming a black body temperature of 0.5 keV, equivalent to an 
upper limit of Zcno < 0.05. 

A third idea for explaining the long tails is that the X-ray 
bursts influence the accretion disk in such a way that the accre- 
tion rate is temporarily increased. A change in accretion rate by 
X-ray bursts has been seen before. For instance, X-ray bursts 
from 4U 1820-303 (Strohmayer & Brown 2002) and 4U 1724- 



307 (Molkov et al. 2000) are so luminous that radiation pressure 
blows away the inner parts of the accretion disk, shutting off 
accretion for a few seconds. Secondly, there are suggestions in 
X-ray bursts from Cen X-4, XTE J1747-214 and 2S 1711-337 
that X-ray bursts act as triggers for switching accretion disks 
from a cold neutral state to a hot ionized state a few days later 
(although that cannot be explained yet in a quantitative manner; 
Kuulkers et al. 2008). Lastly, it appears that immediately be- 
fore and after superbursts the persistent flux from the accretion 
disk behaves differently in the sense that the flux has a some- 
what decreased level for half a day before and an increased level 
for approximately a day after superbursts (e.g., Cornelisse et al. 
2000, 2002; Kuulkers et al. 2002b; Keek et al. 2008). Perhaps 
there is a 4th type of effect from X-ray bursts on disks result- 
ing in long tails. However, 'our' long-tailed X-ray bursts are less 
energetic than the aforementioned bursts that are either super- 
Eddington (i.e., relatively high flux) or superbursts (high flu- 
ence). The X-ray bursts from GS 1826-24 are not particularly lu- 
minous (none show photospheric radius expansion) nor extraor- 
dinarily fluent (Galloway et al. 2006). The same applies to EXO 
0748-676 (only one burst shows photospheric radius expansion; 
Wolff et al. 2005). The majority of the bursts from 4U 1728-34 
do show radius expansion, but have very small fluence due to 
a very short duration (Galloway et al. 2003). Finally, the power 
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Fig. 12. Dynamic depth profile of the specific nuclear energy genera- 
tion rate, as calculated in model 'A3' (Heger et al. 2007). Each level of 
shading indicates a change in the rate by one order of magnitude. The 
bottom of this plot is taken to be the bottom of the reservoir of ashes 
in the model. The hatched regions indicate convective layers (green 
hatched for convection and red cross-hatched for semiconvection; see 
Woosley et al. 2004). The convective region between 10 2 and 10 3 s is 
located at the bottom of the hydrogen left over after the burst and is 
probably due to a Rayleigh-Taylor instability related to composition in- 
version resulting from the burning. It does not have a major influence 
on the structure and evolution of the model. '/JCNO' refers to a variety 
of the CNO cycle where the reaction rates are limited by f} decays; this 
is the 'hot' CNO cycle, '(cf, y)' refers to a capture by heavier isotopes 
such as 12 C and 16 0. A series of bursts is shown; the next burst occurs 
at log(time/s)=4.15. 



density spectra before the burst and during the tail are not signif- 
icantly different (cf., Fig. 6), suggesting no change in accretion 
stream or rate. We believe that an explanation for the long tails 
in terms of a changed accretion environment is less likely. 

5.2. Lack of strong cooling, decreasing black body area 

The lack of cooling could be explained by Compton up- 
scattering of all photons by the hot plasma that surrounds the 
NS and the temperature inferred from the spectrum may in fact 
be representative of the plasma rather than the NS. A problem 
of this explanation is that the inferred emission area drops so 
suddenly in both GS 1826-24 and EXO 0748-676, after 100 and 
60 s respectively. The temperature of the burst at that time is 
still sufficiently high that one should detect large numbers of un- 
scattered photons. A decline of the normalization, if at all, is 
expected to be more gradual. 

Another, in our opinion more likely, explanation is that the 
neutron star may already be fairly hot without the heating by 
thermonuclear flashes. The effective temperature can then never 
drop below the 'quiescent' NS value. There is sufficient persis- 
tent energy production to sustain a NS hot enough to explain our 
measurements, by stable hydrogen burning via the hot CNO cy- 
cle (the model depicted in Fig. 1 1 predicts about lxlO 35 erg s~'), 
pycnonuclear reactions and electron capture processes in the 
crust (the same model assumes 1.6 x 10 34 erg s _1 or 0.15 
MeV/nucleon). Gravitational energy release by the settling of the 
accreted processed matter in the envelope is negligible for these 
accretion rates (Brown & Bildsten 1998). Gravitational energy 
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Fig. 13. Average time profile of 56 X-ray bursts from GS 1826-24 as 
measured with PCU2 on RXTE/PCA in 2-9 keV, now zooming out to 
times further away from the burst. 



release by accretion occurs just outside the NS, is radiated away 
from the NS (King 1995), and is decoupled from the burst emis- 
sion as long as the burst flux is less than the Eddington limit. 
The total energy production rate depends on the mass accretion 
rate and the H and CNO abundance of the donor atmosphere. 
Also, it may partly heat up the core instead of the photosphere. 
Luminosities are expected to reach up to at least a few times 
10 35 erg s _1 . For a canonical 10 km NS radius and ignoring 
the gravitational redshift of a few tens of percents, the Stefan- 
Boltmann law predicts for a luminosity of 2xl0 35 erg s _1 an ef- 
fective temperature of the non-bursting NS of order 0.4 keV. If a 
burst occurs, the photosphere temperature rises to a few ke V and 
the emission is completely dominated by this extra heating, but 
if the temperature becomes comparable to the quiescent value, 
in the tail of the burst, the spectrum is strongly affected by the 
already hot NS. In a standard time-resolved spectroscopic anal- 
ysis, where a pre-burst spectrum is subtracted and the net spec- 
trum is modeled through a black body, one actually subtracts 
one Planck function from another. The resulting function is not 
a Planck function and if, nevertheless, it is modeled as such, the 
physical meaning of inferred emission areas is lost. This effect 
has been extensively studied by van Paradijs & Lewin (1986). 
They find 1) that it is particularly important in the tail of an X- 
ray burst; 2) that the fit results in constant temperatures and de- 
creasing emission areas, as we find in our analyses, and 3) that 
the derived temperature has a tight relationship with the NS tem- 
perature: the measured temperature is about 30% higher than the 
NS temperature outside bursts. Taken at face value, our measure- 
ments imply a NS temperature of 0.7-0.8 keV for both GS 1826- 
24 and EXO 0748-676, ignoring again gravitational redshift and 
deviations from black body radiation that become more signifi- 
cant towards lower temperatures (e.g., Zavlin et al. 1996). This 
is encouragingly close to the simple prediction done above of 0.4 
keV. A similar effect in reversed time order may be happening 
during the rise phase of bursts from GS 1826-24. 

Kuulkers et al. (2002a) investigate this effect in detail for 
bursts from the high-M system GX 17+2. van Paradijs & Lewin 
(1986) suggest to study burst data without subtracting the pre- 
burst spectrum and employing a model that includes components 
for the accretion disk emission and a black body for the thermal 
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emission from the NS. Kuulkers et al. (2002a) follow this sug- 
gestion and model the accretion disk emission by a cutoff power 
law. They find unacceptable values for the goodness of fit and 
dismiss the explanation by an already hot NS and put forward 
the possibility that the decreasing black body radius is connected 
to blanketing effects in the NS atmosphere and comptonization 
of burst photons in the NS atmosphere. Theoretical calculations 
(London et al. 1986; Ebisuzaki 1987; Pavlov et al. 1991; Zavlin 
et al. 1996; Majczyna et al. 2005) show that the color tempera- 
ture is between 1.2 and 1.7 times the effective temperature. The 
fitted black body radius thus decreases by the square of that, to 
maintain the same bolometric flux. This cannot explain the drop 
in radius that we observe which is at least a factor of 5. 

As Kuulkers et al. (2002b) point out, the reason for the un- 
succesful modeling for GX 17+2 is that the accretion disk spec- 
trum also contains a strong black body component due to the 
high M in this LMXB. The disk black body has an only slightly 
higher temperature than that expected of the NS and therefore 
the latter is difficult to distinguish. Also, in GX 17+2 the low- 
energy absorption is high with A^h = 1.9 x 10 22 crrT 2 (Farinelli 
et al. 2007) so that it is tough to find evidence for a NS of tem- 
perature kT « 0.5 keV. Finally, the flux is expected to be low, 
of order 0.1% of the burst peak, so that accumulating a statis- 
tically relevant spectrum is challenging. Our measurements of 
GS 1826-24 do not suffer from these difficulties. Perhaps a de- 
creasing radius and a flattening temperature in a 'standard' burst 
analysis (i.e., with subtraction of the pre-burst spectrum) consti- 
tute the best possible evidence for a hot NS. 

We note that many X-ray bursts in the RXTE catalog 
(Galloway et al. 2008) show similar behavior: temperatures re- 
maining above wl keV and decreasing fitted radii. Since an ex- 
planation by a NS that is already hot without flashes is more 
likely to be applicable to a major portion of the burster popula- 
tion than scattering in a hot circumstellar plasma, this provides 
additional support to that explanation. However, these data are 
vulnerable to lack of low-energy coverage of the PCA, as well 
as from sometimes high absorption columns. As a result, it is 
quite difficult to accurately measure kr below 1 keV. GS 1826- 
24 is one of the few cases where this problem does not exist: we 
have the XMM-Newton data to corroborate the RXTE data and 
a low A^h- 

We furthermore note that the behavior of burst tails in 
UCXBs is expected to be markedly different in the hot-NS sce- 
nario for two reasons. Firstly, UCXBs have much lower H abun- 
dances so that hot CNO burning provides much less heating 
outside flashes. Secondly, in many UCXBs the accretion rate is 
lower so that again the energy production rate is lower outside 
bursts and the NS cooler. This expectation is in line with high- 
sensitivity spectroscopy of some X-ray bursts from UCXBs, for 
example in A 1246-588 (in 't Zand et al. 2008) where the tem- 
perature is seen to decay to 0.5 keV. 

5.3. Long tails in other sources 

Figure 14 shows the time profiles of the bolometric flux for 4U 
1728-34 and EXO 0748-676, together with that for GS 1826-24. 
These are the same data as shown in the top panels of Figs. 4, 
7 and 9, except that one earlier time bin is shown. For conve- 
nience the fluxes have been normalized to the respective peak 
values. This figure shows that there is one property that distin- 
guishes GS 1826-24 from the other two sources: there is a clear 
difference between the prompt and the tail phase. The tails of 4U 
1728-34 and EXO 0748-676 are smooth extensions of the initial 
burst phase. This may be related to the energetic contribution 




Time [sec] 

Fig. 14. The average time profiles of the bolometric flux for GS 1826- 
24 (solid green curve) and 4U 1728-34 (dotted black) and the profile for 
the triple burst from EXO 0748-676 (dashes red; the profile is cut at the 
time the 2nd burst occurs). Fluxes are renormalized with the respective 
peak values. 



of the rp-process being larger in GS 1826-24 (cf, Fig. 12). It is 
probably not related to a smaller amount of inward heating in 4U 
1728-34 and EXO 0748-676. The alternative explanation, deeper 
ignition depths in 4U 1728-34 and EXO 0748-676, is not consis- 
tent with the small burst recurrence times, similar to GS 1826-24 
(Galloway et al. 2008). 

The other documented cases (Swank et al. 1977; Czerny 
et al. 1987; Gotthelf & Kulkarni 1997; Chenevez et al. 2006, see 
Table 2) do seem to have dual time profiles. Their higher tail-to- 
peak flux ratio indicates that the heating to the deeper layers can 
be even more efficient than seen in GS 1826-24. Perhaps there is 
a smaller abundance of heavy isotopes. However, this is difficult 
to assess quantitatively because a quantitative comparison with 
theoretical models has yet to be carried out. 

An exceptional case is the long tail in GX 3+1 (Chenevez 
et al. 2006). It contains much more fluence than the prompt burst, 
perhaps 40 times as much 1 ! Aql X-l (Czerny et al. 1987) and the 
source in M28 (Gotthelf & Kulkarni 1997) may also have large 
fluence ratios but there are data gaps that preclude verification. 
If interpreted as direct NS emission, this cannot be explained 
through cooling of deep layers, particularly since the prompt flu- 
ence is similar to that for other bursts from GX 3+1 that do not 
show a long tail. One difference that distinguishes these other 
cases from the cases discussed in the paper is that they were most 
probably in a different burst regime when the long-tailed bursts 
occurred, namely in a regime without continuous hot-CNO hy- 
drogen burning. 

Several classes of long X-ray bursts, with e-folding decay 
times in excess of roughly 100 s, have been discovered in the 
past decade. This includes superbursts (flashes of 100 m thick 
carbon-rich layers; for a recent list, see Keek & in 't Zand 2008) 
and intermediately long bursts that may result from flashes of 
10 m thick helium layers (in 't Zand et al. 2005; Cumming et al. 
2006; in 't Zand et al. 2007) or otherwise (Chenevez et al. 2007; 
Linares et al. 2008). With the present paper, one realizes that 
even the 'classical' short X-ray bursts can be similarly long in 
some sense. The distinguishing factor is that these bursts initiate 

1 this value depends on where the boundary between prompt burst 
and tail is chosen, but it is always significantly larger than 1 
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with a normal-sized X-ray burst and that, probably, the implied 
thick layer is not heated locally by nuclear reactions but by con- 
duction from a hotter layer on top. 

6. Conclusion 

We have detected an hour-long tail to bursts in GS 1826-24 with 
fluxes and fluences that are two orders of magnitude smaller than 
those of the bursts themselves. We have found similar tails in 
bursts from 4U 1728-34 and EXO 0748-676, although they are 
less distinguished from the prompt burst emission. While detec- 
tion in other bursters is hampered by varying accretion fluxes of 
similar magnitude, there are reports of individual cases of bursts 
with long tails, most notably in GX 3+1 (Chenevez et al. 2006). 
Model calculations show that the tail in GS 1826-24 can be ex- 
plained by delayed cooling of layers that are up to ten times 
deeper in column density than where the flash occurs and that 
were heated up through inward conduction of the flash heat. 
Possibly tails in other sources can be similarly explained. Further 
model calculations are needed, where dependencies of donor 
composition, ignition depth and accretion rate are taken into ac- 
count. Comparing such calculations to different kinds of bursts 
and bursters may yield constraints on the details of conduction 
in the neutron star envelope. 

A characteristic of the tails, that at first hand is unexpected in 
this scenario, is the small amount of cooling. Rather than being 
due to up-scattering of the burst photons by a hotter optically- 
thick plasma, we believe that the most likely explanation is that 
the NS is already hot without flashes. The temperature that one 
measures in the burst tail is then representative for the non- 
bursting NS. This scenario also provides a more natural expla- 
nation for the decreasing black body normalization in the tails. 
As discussed by van Paradijs & Lewin (1986), this presents an 
interesting opportunity to study NS temperatures as a function 
of accretion rate. However, this may be a cumbersome. Apart 
from good low-energy coverage, one would need to work on a 
single-burst basis, also for the accurate modeling of the persis- 
tent spectrum. Better data would be needed than presented here. 
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